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On the broadening of emission lines in active galactic nuclei 
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Abstract. Recent works, both theoretical and observational, have suggested that turbulence could play a non- 
negligible role in the broadening of emission lines in active galactic nuclei. The purpose of this note is to show 
how shock wave-turbulence interaction, under unsteady regime, can affect the broadening of emission lines. 
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1. Introduction 

The Seyfert galaxies are a class of active galactic nuclei 
(AGN) characterized by nuclear spectra with strong emis- 
sion lines by highly ionized atoms. These lines are of two 
types: narrow, corresponding to low density plasma with 
bulk velocities of the order of 100 km/s; broad, corre- 
sponding to high density plasma and bulk velocities of the 
order of 10 4 km/s. The main physical mechanism responsi- 
ble for the ionization is the non-stellar photon flux coming 
from the AGN itself (for reviews, see Davidson & Netzer 
1979, Krolik 1999, Fabian et al. 2000, Sulentic et al. 2000). 
However, also shocks can provide the physical conditions 
necessary to the ionization of the atoms. This occurs, for 
example, when the postshock temperature reaches values 
higher than 1 eV, and therefore the shocked gas emits a 
continuum in the EUV-soft X-ray band, which is quickly 
absorbed by nearby gas and its energy converted into emis- 
sion lines. See also, however, the work of Perry & Dyson 
(1985), Collin-Souffrin et al. (1988), the former recently 
updated and enlarged by Fromerth & Melia (2001). 

The main objection to shock formation of line emission 
is that the wide range of observed ionization stages is not 
fully compatible with a thermal emission from post-shock 
regions. The temperature required to produce the highest 
stages are far from those obtained in shocks. In addition, 
according to Krolik (1999), we have no direct evidence 
of the existence of shocks, but other authors suggested 
that shocks can be important in the morphology, kine- 
matics, and excitation of emission lines in Seyfert galaxies 
(Colbert et al. 1998, Allen et al. 1999). Indeed, the emis- 
sion line ratios do not demonstrate the excitation mech- 
anism (shocks or otherwise). Moreover, the velocity field 
and line width are consistent with shocks and with en- 



ergetics determined from radio and X-ray observations 
(Allen et al. 1999). 

Turbulence has been invoked to explain the line broad- 
ening (Horne 1995, Pariev & Bromley 1998), but some au- 
thors object that the turbulent velocity field in an AGN 
disk is not so strong to affect the line profile (Fabian et 
al. 2000). Bottorff et al. (2000) have found that extra 
thermal, non-dissipative microturbulcncc somewhere be- 
tween and 10 3 km/s, is consistent with observations of a 
small sample of quasars with high signal to noise spectra 
(Baldwin et al. 1996). In addition, Bottorff & Ferland (pri- 
vate communication; to be published on ApJ, April 1, 2002 
issue) find that ss 200 km/s microturbulence with dissi- 
pation is able to explain typical emission line ratios found 
in quasars and other AGN. Observations poorly constrain 
the range of turbulence in AGN, but even the possibility 
of its presence has important ramifications in the analysis 
of AGN emission lines. It is therefore useful to investi- 
gate mechanisms that generate and sustain (possibly high 
levels of) turbulence. 

The purpose of this note is to discuss how the inter- 
play between shock waves and turbulence can enhance the 
turbulent velocity field. These physical issues have been 
already applied to the line broadening in pulsating stars 
(Gillct et al. 1998) and we try to apply these concepts to 
active galactic nuclei. 

2. Fully developed turbulence 

Before going on, let us to remember some concepts in 
fully developed turbulence. We refer to Landau & Lifshitz 
(1987). Turbulent flow is characterized by random varia- 
tions in space and time in the velocity field and develops 
at very high Reynolds number: 
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where I is the characteristic dimension of the flow, V is the 
velocity, p the mass density, and \x the viscosity, which is 
equal to about KT 17 T 5 / 2 g • cm s for a fully ionized 
unmagnetized hydrogen plasma (see Choudhuri 1998). For 
a typical electron temperature of 10 4 K, we have p = 1CP 7 
g • cm _1 s _1 . 

By considering a typical flow speed of V2 km/s, where 
V2 = V/100 km/s, and a mean density of 1.67 ■ 10~ 15 
g/cm 3 (i.e. an electron volume density of 10 9 cm~ 3 for a 
fully ionized hydrogen plasma), we obtain from Eq. (Q) 
that Re w 0.167 • I ■ V2, where I is in cm. Therefore, when 
the characteristic dimension of the flow is of some hun- 
dreds of metres we have already a turbulent flow (for a 
speed of the order of 100 km/s). On astrophysical scales, 
we can consider that the characteristic length is of the or- 
der of astronomical units (AU), so that the flow is fully 
turbulent, even though in presence of tiny values of vis- 
cosity. 

In a turbulent regime, we can divide the velocity V 
into a mean value component v, obtained from an aver- 
age over a long time, plus a fluctuating part u that when 
averaged is equal to zero: 

V = v + u (2) 

The order of magnitude of the distance over which the 
speed varies appreciably (i.e., the dimension of the largest 
eddies) has the same order of the dimension I of the region 
where there is the flow. Therefore, the fluctuation of the 
speed is of the same order of the variation of the mean 
speed over the scale I. However, under normal conditions, 
the speed fluctuations cannot exceed the local sound speed 
c s (cf. Balbus & Hawley 1998). With respect to time, the 
velocity does not vary appreciably over a timescale smaller 
than about I /v. 

Since the viscosity decreases as the Reynolds number 
increases, in larger eddies there is no or negligible vis- 
cosity, and hence no energy dissipation. The energy flows 
from larger to smaller eddies, where it is transformed into 
heat and dissipated. To maintain a steady state flow, an 
external source of energy must be provided to supply en- 
ergy to large eddies. 

In the plasma of an accretion disk around a supermas- 
sive black hole, in addition to the fluid viscosity, other 
effects take place, because of the nature of the plasma. 
In this case, the electromagnetic forces play an important 
role in transferring energy and momentum. The change in 
state, particularly during the passage of a shock, derives 
from the collective interactions between charged particles 
and electromagnetic fields. The fields can be mainly of 
two types: constant in time (i.e. produced by space charge 
separation, currents, or other external sources) or fluctu- 
ating in time (i.e. produced by plasma instabilities). The 
first case is usually referred as laminar, while the second 
as turbulent. For more details see, for example, Tidman 
& Krall (1971). 



3. Shock waves-turbulence interaction 

Despite its importance in several application of science 
and engineering, the interaction of shock waves with tur- 
bulence is still largely unknown. The first problem is to 
obtain reliable experimental data at high Mach numbers 
[M > 5, hypersonic flow), but also experiments at mod- 
erate Mach numbers suffer with a strong dependence on 
the measurement instruments and therefore cannot de- 
scribe with reasonable reliability the complex interaction 
of shock waves with turbulence. 

Recently, the availability of computer power made it 
easier to set up numerical models, but problems remain, 
because, when dealing with turbulence, it is necessary to 
introduce closure conditions, which in turn are dependent 
on experimental data. For a review, see Andreopulos et 
al. (2000), Lele (1994), Adamson & Messiter (1980), and 
references therein. 

However, there is a general consensus that the most 
important feature of the interaction between shock waves 
and turbulence under unsteady conditions is the amplifica- 
tion of speed fluctuations and a strong change in the scale 
lengths, even though it is not well understood how this 
amplification occurs (Andreopulos et al. 2000). According 
to actual studies, the amplification depends on the shock 
strength, the state of the turbulence, and its level of com- 
pressibility. With respect to the amplification of the tur- 
bulent kinetic energy, Rotman (1991) reports an increase 
of a factor 2 — 2.15, while for Jacquin et al. (1993) the 
amplification depends on the density ratio: 

k 2 {t) 2 + p 2 ft)/p 2 (0) 

k 2 (0) 3 [ ' 

where k(t) and p(t) are the turbulent kinetic energy and 
the gas density at the time t, respectively. In this case, 
the amplification of the kinetic energy can be up to a/6 
for monatomic gases (see the next Section). 

It is worth remembering, that at high Mach number 
and for steady state motion, the effect of compressibility 
suppresses the turbulence and therefore there is no ampli- 
fication. Changes in the Mach number (unsteady motion) 
generate distortions in shock waves, allowing the interac- 
tion with turbulence. The key point is therefore that the 
Mach number must not be constant. According to Jacquin 
et al. (1993), these conditions represent the "pressure re- 
leased" regime and are the upper limit of the amplification. 
In this case, the turbulent Mach number can be higher 
than 1. 

4. Post-shock conditions 

For the purposes of the present note, we rearrange the 
Rankine-Hugoniot relations in the limit of a strong shock, 
that is with Mach number M >> 1 (for reviews see, for ex- 
ample, Zel'dovich & Raizer 1967, Landau & Lifshitz 1987, 
Frank et al. 1992, Krolik 1999). The post-shock pressure 
P 2 is: 

p 2 n(l-^)(l + a) Pl V? (4) 

92 
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where a is the ionization factor (a = 1 full ionization), pi 
is the post-shock density, p\ and V\ are the undisturbed 
density and speed of the plasma. Here, and in the following 
we use the subscripts 1 and 2 to indicate the undisturbed 
and the post-shock values, respectively. 

The density ratio across the shock is given once only 
the value of the specific heat ratio 7 is known: 



P2_ 

Pi 



7 + 1 



7' 



1 



(5) 



From the equation of continuity across a shock, we can 
get easily that the ratio of velocities is a function of 7 only: 



v 2 



P2 
Pi 



7 + 1 

~ T 



(6) 



In studies of AGN, velocities widths of typical bright 
BLR emission lines are non relativistic, therefore 7 is gen- 
erally considered to be 5/3, the value of monatomic gas 
or metal vapors. Therefore, by considering 7 = 5/3 and 
Eq. (H), we obtain: 



P 2 « 0.75(1 + a) Pl V? 



(J) 



To calculate the post-shock temperature, we need of 
the equation of state: 



P = (1 + a)pRT 



(8) 



where R is the gas constant for the considered plasma. 
Therefore, by taking into account the Eqs. (|5|), ([?]), and 
(||), we can arrange the numerical value of the post-shock 
temperature as a function of the undisturbed speed: 



To 



2(7-1) 



(l + a)i?( 7 + l) : 



rn 2 



(9) 



Please note that V\ is given by Eq. , by considering 
the mean value. We could consider also the peak value, 
given by the sum of the average plus the fluctuating tur- 
bulent value. 

5. The broadening of emission lines 

There are several models to explain the broadening of 
emission lines by shocks (Perry & Dyson 1985, Collin- 
Souffrin et al. 1988, Fromcrth & Melia 2001). Recently, 
turbulent motion has been invoked to explain part of the 
broadening (Bottorf et al. 2000), even for Fe Ka fluores- 
cence line (Pairev & Bromley 1998). Herein we apply tur- 
bulence amplification by shocks to investigate where in 
the range to 10 3 km/s (Bottorff et al. 2000) turbulence 
(if it exists in AGN) may lie. 

As the electron temperature of the environment 
plasma was measured at 10 4 K, if we assume that pro- 
tons have the same temperature of electrons, this implies 
a maximum fluctuation speed of some tens of km/s (we 
have considered that under normal conditions the turbu- 
lence is limited by the local speed of sound; cf. Balbus & 
Hawley 1998). 



On the other hand, we can assume that the electron 
temperature is not representative of the mean temperature 
of the gas. If we consider a shock due - for example - to 
a wind with speed, say 500 — 1000 km/s, the post-shock 
temperature, with reference to Eq. (||), is about 1.1 — 4.6 • 
10 7 K (assuming a — 0.5 and a hydrogen plasma). This 
is the ion temperature, while the electron value is scaled 
by a factor m e /m, « 5.4 • 10~ 4 (T e w 0.6 - 2.5 • 10 4 K, a 
value comparable to the observations). Indeed, it is known 
that the basic feature of shock waves in plasma is the 
slow energy exchange between ions and electrons, and the 
high electron mobility. The equilibrium temperature will 
be reached slowly after a time t eq : 



dTe 

dt 



T-T e 

tea 



(10) 



because the energy exchange occurs by Coulomb scatter- 
ing and there is large difference of mass between ions and 
electrons. These relaxation processes determine the thick- 
ness of the shock front (Zel'dovich & Raizer 1967). 

From the post-shock temperature, we can calculate 
the most probable thermal speed, assuming a Maxwcllian 
distribution. Indeed, even though the electrons and ions 
do not exchange energy among them, they reach indepen- 
dently and quickly the Maxwellian distribution (generally 
after one collision). So, the most probable speed is: 



Vth 



2kT 



(11) 



and the value of the turbulent speed, that is limited by 
the local sound speed (cf. Balbus & Hawley 1998): 



u = c. = 



^kT 



(12) 



If the shock interacts with the turbulence under un- 
steady conditions, we have the amplification of the turbu- 
lent speed. We know from Sect. 3, the amplification of the 
kinetic energy, from which we can calculate the amplified 
speed, under the simplified assumption that the distortion 
affects only the speed and not the density. We introduce 
the speed amplification factor (3 that is simply the square 
root of the factor for the kinetic energy. According to the 
theory by Rotman (1991), (3 = while for Jaquin et al. 
(1993) it depends on the density ratio across the shock: 



(3 = 



2 + pi/ pi 



(13) 



and can reach the value of 1.57 for a monatomic gas. 

Therefore, the maximum value of amplified turbulent 
speed is: 



(14) 



Substituting Eq. (g) in Eqs. (0), ((T|), and (0), we 
obtain that V th « 0.87Vi, u « O.8V1, and 1.41 Vy < u a < 
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Fig. 1. Evaluation of post-shock thermal (A), turbulent (B), and amplified turbulent speed (C, according to Jaquin 
et al. (f993); D, according to Rotman (199f )) as a functions of the incoming shock speed. 



1.57Vi, according to the theory by Rotman or Jaquin 
et al., respectively. Some order-of-magnitude calculations 
are shown in Fig. |l|. 

The maximum value of the Doppler broadening Ap/uq 
of the emission lines due to the combined action of thermal 
and amplified turbulent speed is (see Rybicki & Lightman 
1979): 



Av 1 kT, n no . 
— = ~J— 2 + /P-y) 
fo cynii 



(15) 



where c is the speed of light. An evaluation of the in- 
fluence of the amplification of turbulence on the Doppler 
broadening is shown in Fig. || for hydrogen plasma. 

We can also do the reverse calculation, i.e. starting 
from post-shock temperature we can infer the speed. For 
example, we know that bright emission lines stop forming 
efficiently when T w 10 5 K. For a fully ionized (a = 1) hy- 
drogen plasma, Eq. (||) gives the pre-shock speed, which 
has the values of about 94 km/s. For the amplified tur- 
bulent speed, we obtain a value between 133 and 148 
km/s, according to Rotman or Jaquin's theory, respec- 
tively. These values can be lower for lower ionization: if 
a = 0.5, we have that the mean value of the speed is be- 
tween 115 and 128 km/s. Note that values change if we 
deal with heavy ions. 

Several authors have computed models for spectra 
from shocks (e.g. Cox 1972, Daltabuit et al. 1978, Shull & 
McKee 1979, Binette et al. 1985) under the assumptions 
of steady state and constant pressure in the post-shock re- 
gion. These conditions do not hold anymore in unsteady 
state and it would be necessary to set up a numerical 
model to evaluate the impact of the amplification of the 
turbulence on the emitted spectrum. However, we can now 



make some inferences: radiation generated into the shock 
acts to smooth the amplification, so as to establish again 
a sort of equilibrium. We know that hydrogen line intensi- 
ties are very nearly proportional to pv 3 (Cox 1972). In this 
note, we assumed, for sake of simplicity, that turbulence 
affects only the speed and not the density, so we expect 
an increase of the hydrogen line intensities due to the am- 
plification of the turbulent speed. Moreover, because the 
post-shock temperature can be high enough (w 10 7 K), 
we expect also an increase of X-rays emission by thermal 
bremsstrahlung (if the shocked plasma is thin; otherwise 
a blackbody temperature is detected). 

It is worth noting, that these are order-of-magnitude 
calculations made with the speed modulus, but we can ob- 
serve and measure the radial component of these vectors. 
So we would introduce in calculations also the line-of- 
sight angle. 



6. Final remarks 

We emphasize that the purpose of this research note is to 
settle basic concepts about the turbulence amplification 
in AGN and some of the consequences, mainly related to 
the broadening of emission lines. Although, in this case 
it is an additional term only, the impact of this physical 
mechanism on the observations is still to be determined, 
particularly with reference to LINER nuclei (cf. for ex- 
ample, Dopita et al. 1997, Allen et al. 1999, Evans et al. 
1999). 

The theory suffers from the lack of understanding of 
how amplification occurs - as underlined by several other 
authors - and therefore, for this application, observations 
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Fig. 2. Evaluation of broadening due to thermal speed only, additioned with turbulence, and amplified turbulence 
(according to Jaquin et al. 1993, i.e. the upper limit). Values are calculated for an incoming speed of 1000 km/s. 



of some AGN in soft X-rays or EUV are required to better 
assess the problem. 
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